Distinct order of Gd 4/ and Fe 3d moments coexisting in GdFe4Al8 
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Single crystals of flux-grown tetragonal GdFe4Al8 were characterized by thermodynamic, trans- 
port, and x-ray resonant magnetic scattering measurements. In addition to antiferromagnetic order 
at Tjv « 155 K, two low-temperature transitions at Ti « 21 K and r2 « 27 K were identified. The Fe 
moments order at Tjv with an incommensurate propagation vector (r, r, 0) with r varying between 
0.06 and 0.14 as a function of temperature, and maintain this order over the entire T < Tjv range. 
The Gd 4/ moments order below T2 with a ferromagnetic component mainly out of plane. Below 
Ti, the ferromagnetic components are confined to the crystallographic plane. Remarkably, at low 
temperatures the Fe moments maintain the same modulation as at high temperatures, but the Gd 
4/ moments apparently do not follow this modulation. The magnetic phase diagrams for fields 
applied in [110] and [001] direction are presented and possible magnetic structures are discussed. 

PACS numbers; 75.25.-|-z, 75.30.-m 



I. INTRODUCTION 

Understanding the interplay of rare-earth local- 
moment magnetism and 3d transition metal itinerant 
magnetism is of fundamental physical interest and may 
help in the design of more efficient permanent magnets. 
The tetragonal i?Fe4Al8 {R =rare earth) compounds 
(Fig. nj are well suited for studying this interplay be- 
cause of simple symmetry conditions and because the in- 
teraction between the two magnetic sublattices is rather 
weak: Fe moments appear to order without a corre- 
sponding R 4/ order. Consequently, numerous studies 
have been performed on i?Fe4Al8 (R148) over the last 

30 years and interest in these compounds has remained 
^^jg^ 2.3.4.5.6.7.8.9.10.11 

Neutron scattering studies performed on R148 with 
various R indicate that Fe moments, between which 
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FIG. 1: (Color online) The crystal structure of i?Fe4Al8. Thin 
lines outline the tetragonal (I4/mmm) unit cell, thick lines 
denote the Fe "cage" around the rare earth R. 



the couphng is strongest, order between 100 and 200 K, 
generally in cycloidal structures with propagation vec- 
tors parallel to [110] and moments confined to the a—b 
plane i^i^iiSiii For magnetic rare earths R, ordering of the 
R moments, with much weaker coupling and at much 
lower temperature T, has been reported. However, the 
critical temperature associated with this ordering is very 
poorly defined. For example in Dyl48 at ~50K the 
magnetic dc susceptibility rises faster than expected for 
Curie- Weiss behavior, but there are no sharp features!^ 
Furthermore, the rare earth moments were found to fol- 
low the Fe-moment modulation^iifl This seems to imply 
that the Fe— i? moment interaction is stronger than the 
R—R moment interaction. 

Of all rare earth elements, Gd has two highly inter- 
esting specific features associated with its high-spin and 
zero-angular- momentum 4/ state: 1) high spin results 
in strong magnetic coupling between the localized 4/ 
and the conduction electrons, implying large magnetic 
interactions and the largest de Gennes factor of all rare 
earths; 2) zero angular momentum implies a spherical 
4/ charge-cloud and no magneto-crystalline anisotropy 
(MCA) resulting from crystal-electric-field effects — the 
direct interplay of magnetic interactions can be studied 
without crystal-electric-field effects. In view of this, it 
is somewhat surprising that GdFe4Al8 (Gdl48) has been 
much less studied^'^ than R148 with other rare earths. 
In particular, apart from Mossbauer spcctroscopyjSi there 
have been no microscopic studies of Gdl48. 

Here, we report on the flux-growth of single crystals 
of the GdFe4Alg phase, and on the characterization of 
the crystals by magnetization, electrical transport, spe- 
cific heat, and x-ray diffraction. We also report on first 
synchrotron x-ray resonant magnetic scattering (XRMS) 
data. We provide evidence for two consecutive phase 
transitions at low temperature T in addition to the Neel 
transition at 155 K, and show that the two, low-T, tran- 
sitions are associated with the ordering of the Gd 4/ 
moments, resulting in complex magnetic structures in- 
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volving both ferro- and antiferromagnetic components. 
In contrast to other R148 compounds, the 4/-moment 
ordering is associated with a propagation vector that is 
distinct from that associated with the Fe moment order, 
which at low T has the same modulation as in the vicinity 
of the Neel temperature. We also present the magnetic 
phase diagram for two field directions and discuss the 
implications of changes in the Fe stoichiometry. 

The paper is organized into seven sections. In Sec. 
lU we describe the flux growth of the single crystalline 
samples used in the study, and the experimental proce- 
dures. In Sec, mil zero and low field electrical transport, 
magnetization, and specific-heat data are presented, and 
we evaluate the influence of the iron stoichiometry on 
the physical properties. In Sec. II VI we present electrical 
transport and magnetization data measured in fields for 
two field directions, and in Sec. the results of a first 
XRMS experiment are described. Finally, we discuss the 
low-temperature magnetic phase diagram and possible 
magnetic structures in Sec. IVII before summarizing our 
main conclusions in Sec. IVIII 



II. EXPERIMENTAL 

Whereas all previous single crystals of _RFe4Alg ma- 
terial were obtained using the Czochralsky-method, we 
have grown single crystals of Gdl48 (and other R148) 
with a self-flux methodiiSiiii^ One of the problems with 
studies on single crystals of R148 is a width of formation 
often observed, involving some Fe atoms occupying nom- 
inal Al positions and vice versa; this can lead to striking 
differences in the magnetic structure and phase diagram 
as com par ed to the stoichiometric material (see, e.g., 
Refs. I9II15I). A flux-growth procedure may allow better 
control of the stoichiometry of the crystals by varying the 
starting composition. Unfortunately, information about 
the ternary phase diagram i?-Fe-Al is very limited. We 
used differential thermal analysis (DTA; see Ref . IT^ for a 
review of the use of DTA in the flux growth of crystals) to 
establish i) that GdFe4Al8 is congruently melting and ii) 
selected solidus temperatures in the ternary around the 
GdFe4Al8 stoichiometry. For the crystals used in this 
study, we selected starting compositions between Gdl48 
and the also congruently melting Fe2Al5 (Gd3Fei6Al34 
and Gd2Fei4Al3i), which should not lead to large imbal- 
ances in the Fe to Al ratio. 

Elements in ratios corresponding to the starting com- 
position were first arc-melted together. The resulting 
ingot was placed in an alumina crucible, which was 
wrapped in Ta foil (to prevent any residual oxygen in 
the Ar from reacting with the sample) and placed in a 
vertical tube furnace in flowing Ar. Crystals were grown 
by heating to 1475°C and then slowly (2°C/h) cooling to 
1180°C, at which point the furnace was turned off. The 
flux was removed from the crystals in a second step by 
heating to 1200°C, keeping the temperature for 30min 
and then decanting, following procedures described in 



Refs. ll2ll3ll1ll 

For both starting compositions, we obtained crystals of 
the Gdl48 phase, as identified via powder x-ray diffrac- 
tion. Crystals typically grow prismatically with (110) 
facets and the long direction parallel to [001], as de- 
termined with Laue scattering. We obtained crystals 
of dimensions up to about 10 x 2 x 2mm'^. We deter- 
mined, via single crystal x-ray-diffraction-structure re- 
finement and electron-microprobe analysis (employing 
single crystals of GdFe2 as standards), the crystals to 
be slightly iron deficient, but detected no Al on Fe 
sites or vice versa (measured compositions were between 
GdFe3.88(5)Al8 and GdFe3.g6(i) Al8; deviations of the Al 
stoichiometry from 8 were always less than 2 standard 
deviations and are not listed). 

We note that the x-ray-diffraction-structure refinement 
cannot distinguish between a Fe deficiency due to vacan- 
cies and a (larger) Fe deficiency due to Al on Fe sites, 
both scenarios can give an equal electron density. How- 
ever, assuming an Al/Fe mixture to retain the full oc- 
cupancy of the Fe site leads (for one of the investigated 
crystals) to a GdFe3 92(1) Al8.o8(i) composition and, thus, 
to an excess of Al. Since the electron-microprobe analysis 
results indicated no significant variations in the 1/8 ratio 
of Gd to Al, we did the structure refinements under the 
assumption that the electron density on the Fe sites is 
caused by vacancies, although some Al/Fe substitution 
could never be excluded. Whatever scenario is chosen, 
none of our conclusions drawn in this article are affected 
by the issue of whether or not some Al atoms are present 
on the Fe site, and for definiteness subsequently we will 
assume that the vacancy scenario is the correct one. The 
composition of our samples can thus be described with 
the empirical formula GdFe4_5Al8, with 5 « 0.04 — 0.12. 
Surface scans in the electron-microprobe analysis gave no 
indications of compositional variations within the same 
crystal. 

Since the exact iron stoichiometry varies slightly from 
crystal to crystal (even within the same growth batch) 
and this was found to have a significant influence on 
the magnetic properties (see Sec. IIIII below 1 most of the 
measurements presented below were performed on the 
same crystal, sample I, which has a refined composition 
GdFe3.96(i)Al8 (a = 8.7699(9) A, c = 5.0440(6) A). A bar 
cut from the crystal by a wire saw was connected with 
contacts for 4-point electrical-transport measurements 
with the current density j||[110] (sample la). Other sam- 
ples were connected with contacts for electrical transport 
measurements with either j||[110] or j||[001]. On a sec- 
ond bar (sample lb, which had a mass of 2.22 mg) cut 
from crystal I we measured the longitudinal magnetiza- 
tion and the zero-field specific heat. In order to inves- 
tigate the relation between compositional variations and 
variations in the transition temperatures, we carried out 
additional magnetization measurements on several other 
crystals, the composition of which was determined by x- 
ray diffraction and electron microprobe. The additional 
crystals are labelled II to VII in the order of their ap- 
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pearance in the text. The above measurements were per- 
formed with commercial (Quantum Design) laboratory 
equipment. For the high-field measurements, the sample 
orientation was fixed with a two-component glue (X60 
from Wagetechnik GmbH, Darmstadt, Germany). 

Sample lb was also used in the XRMS experiment, per- 
formed on the 6ID-B undulator bcamlinc in the MUCAT 
sector at the Advanced Photon Source, Argonne National 
Laboratory. The incident energy of x-rays was tuned to 
the Gd Lii edge {E = 7.934keV), using a liquid nitro- 
gen cooled, double crystal Si (111) monochromator and a 
bent mirror. Sample lb was mounted on a copper rod on 
the cold finger of a closed cycle displex refrigerator, such 
that a natural (110) facet was exposed to the x-ray beam. 
Thermal transfer was enhanced by embedding the sam- 
ple in copper paste. The sample, oriented so that the 
scattering plane of the experiment was coincident with 
the a — b plane, was encapsulated in a Be dome with He 
exchange gas to further enhance thermal transfer. The 
incident radiation was linearly polarized perpendicular 
to the scattering plane [a polarized). In this geometry, 
only the component of the magnetic moment that is in 
the scattering plane contributes to the resonant magnetic 
scattering arising from electric dipole transitions (El) 
from 2p to 5d states. Furthermore, the dipole resonant 
magnetic scattering rotates the linear polarization into 
the scattering plane (tt polarization). In contrast, charge 
scattering does not change the polarization of the inci- 
dent beam {a ^ a scattering). Pyrolytic graphite PG 
(006) was used as a polarization analyzer, selecting pri- 
marily TT polarized radiation. For E = 7.934 keV, the 
polarization analyzer used reduces the detected inten- 
sity resulting from a ^ a charge scattering by about 
99.9%, whereas the cr — * tt resonant magnetic scattering 
is passed with little loss. Thus, the polarization analysis 
suppresses the charge scattering relative to the magnetic 
scattering signal. 



III. LOW FIELD MEASUREMENTS AND 
INFLUENCE OF FE CONTENT 

In this section, we will first present zero-field transport 
and low field magnetization measurements on our main 
crystal (samples la/Ib). Then, we will briefiy discuss 
sample-to-sample differences in the characteristic tem- 
peratures, which we think are related to Fe deficiency. Fi- 
nally, we will present specificheat measurements on sam- 
ple lb indicating that the Gd 4/ moments are ordered 
only at low temperatures. The main implications of the 
measurements will be discussed later in Sec. I VII 

In Fig. 121 we show, for sample la, the T dependence 
of the resistivity p (j||[110]). Since the sample is rather 
small, the estimated uncertainty in the absolute value of 
the resistivity is about 15%. The antiferromagnetic or- 
dering at Tn — 155 K is visible as a sharp kink in p(T)m^ 
The inverse dc magnetic susceptibility (M/H) of 
sample lb, cut from the same crystal as sample la, is 
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FIG. 2: (Color online) Zero-field resistivity p (along [110], 
sample la) and inverse dc susceptibility vs temperature 
T (sample lb). 



also presented in Fig. |21 The dashed and dotted curves 
shown are determined from the magnetization M(T) in a 
field oi p,o H — O.n applied parallel to [110] and parallel 
to [001]. Below this field M vs H is finear for T > 40K. 
Since x is slightly anisotropic, we also calculated the 
polycrystaUine average Xpoiy = (2xiio + Xooi)/3, shown 
as squares. In for both directions the antiferromag- 
netic ordering is manifest by a change in slope, but it 
is much less sharp than the corresponding signature in 
the resistivity. A weak signature of the Neel transition 
may be expected if, as we will discuss below, only the 
Fe moments order at T/v. In agreement with measure- 
ments on polycrystaUine samples,^ x follows Curie- Weiss 
behavior both above and, over a limited T range, below 
T/v. We found (on the calculated polycrystaUine average) 
above T/v a Weiss temperature of —165 K and an effective 
moment of 11.4 /f.u.. Assuming contributions by the 
magnetic atoms that are additive in the Curie constant 
and a contribution by Gd of 7.9 /x_b /Gd, the free moment 
per Fe atom in the paramagnetic state is 4.1 /is /Fe. Be- 
low Tpf (range 70 — 115 K) we found a Weiss temperature 
of 17 K and an effective moment of 7.4 pB /f.u.. The lat- 
ter is closer to the Gd free-ion value {7.9 pB /f-u.) than 
the one reported by Buschow et al? {6.2 ps /f.u.). 

The resistivity curve shows additional features just be- 
low 30 K, magnified in Fig. |3i) (thick line): The broad 
peak in p{T) at ~28K with a drop at ~26.5K (T2) 
indicates a phase transition. Furthermore, with decreas- 
ing T, at ■~21K (Ti) there is a sudden increase in p by 
^ 10%. This feature is hysteretic in temperature, sug- 
gesting that the feature is associated with a first-order 
transitionii^ Note that the hysteresis visible at T2 is op- 
posite to what is expected for a first-order transition (i.e., 
superheating/supercooling). It is most likely a remnant 
of the hysteresis associated with the Ti transition. At 
Ti, the form of the hysteresis is consistent with super- 
heating/supercooling. 

Also shown in Fig. |3^) are low-T zero-field resistivity 
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FIG. 3: (Color online) Low-temperature normalized resis- 
tance (panel a) and derivative of high temperature resistance 
(panel b) vs temperature T of various crystals: la (thick line) 
and 11 (dotted line) for j||[110], III (open circles) and IV (full 
circles) for i||[OQl]. 



curves for three additional samples, one also with j\\ [110] 
and two with j||[001]. Whereas the feature of a sudden 
increase of p with decreasing T is visible in all curves, 
the other feature at slightly higher temperatures is only 
identifiable when j|| [110]. Based on the results from sam- 
ple la only, it may be tempting to associate the increase 
in p (with decreasing T) with the opening of a super- 
zone gap^ However, the presence of this increase with 
the same order of magnitude for both j|| [110] and j\\ [001] 
makes a superzone gap scenario less likely. 

It is also clear from Fig. that there are significant 
variations of the temperatures Ti corresponding to this 
feature. For comparison. Fig. (Jja) shows the resistivity 
derivative dp/dT for the same samples at higher T, in 
order to make the Neel transition more visible, which 
manifests itself differently for the two current density di- 
rections. The curves indicate a sharp, well defined Tjv for 
each sample and comparing Figs.|3tt) andlSb) suggests a 
correlation between higher Ti and lower Tjy. 

Figure 21 displays low field (iJ||[001]) magnetization 
measurements performed on three samples with differ- 
ent Fe stoichiometry as determined by x-ray diffraction 
and electron microprobe: samples lb (GdFcs 95(1) Alg, 
thick line), V (GdFcs 94(g) Alg, full squares), and VI 
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FIG. 4: (Golor online) Low field (i/||[001]) Magnetization of 
various crystals with Fe stoichiometry (as described in Sec. 
rm determined by electron microprobe or x-ray-diffraction- 
structure refinement. 



(GdFe3.88(5) Alg, dashed line). The data were taken for 
both decreasing and increasing T (field cooled), and 
normalizedSS to the maximum in magnetization. A com- 
parison with Fig. indicates that the rise in M with 
T is associated with Ti, the subsequent decrease with 
T2. A decrease in the Fe content seems to systematically 
shift Ti and T2 to higher temperatures. We note that for 
crystal I M [T) (Fig. ^ indicates a higher T2 than p{T) 
[Fig. Oi)]. In principle, this might be due to a different 
Fe content of the two pieces (la and lb) cut from the 
crystal and used for these measurements. However, we 
do not believe that this is the case, because the specific- 
heat peak associated with T2 (see below; measured on 
lb) corresponds well to T2 deduced from p{T) on sample 
la. 

Combining the above data it seems that decreasing the 
Fe content results in a decrease of Tjv and an increase of 
Ti and T2. In order to confirm this, a study involving 
samples that are much more Fe deficient, would be highly 
desirable. 

Figure^lshows the measured specific heat Cp of sample 
lb (closed squares) after subtraction of the addenda con- 
tribution (the specific heat of sample platform and grease 
measured before mounting the sample). Also shown is 
the mass-scaled specific heat of YFe4Al8 from Ref.'2l'(full 
line), and the difference between the two specific heats 
(open circles) . Since the necessary mass-scaling is small, 
the remaining specific heat is close to the magnetic con- 
tribution associated with the Gd sublattice. The broad, 
asymmetric, peak at ~26.5 K corresponds to the drop in p 
at T2 [determined on sample la cut from the same crystal 
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FIG. 5: (Color online) Specific heat Cp measured on sample 
lb in zero field (■) and Gd 4/ moment contribution to the 
specific heat (o), obtained by subtracting the specific heat of 
YFe4Al8 (full line, from Ref.E3). Inset: Entropy of the Gd 4/ 
moments obtained by integrating the 4/ moment contribution 
to Cp/T. The two lines shown were calculated using two 
extrapolations of the specific heat to T = 0, indicated by 
dashed lines in the main panel (see text for details of the 
extrapolations) . The dotted line in the inset indicates the full 
Gd 4/ entropy. 



as sample lb, see Fig. Ot) dashed curve]. The broadness 
of the specific-heat peak may suggest a degree of inho- 
mogeneity of the Fe distribution in the sample, although 
the electron-microprobe analysis provided no indications 
of inhomogeneities within crystals. 

No feature in the specific heat is visible around Ti, 
the small discontinuity at 20 K indicated by the data was 
found to be an artifact due to the change of the heat pulse 
intensity applied by the system. We checked the specific- 
heat raw data, and can exclude any latent heat restricted 
to a temperature region smaller than the spacing of mea- 
surement points. We conclude that the latent heat of the 
Ti transition has to be small, i.e., the phases above and 
below Ti have similar entropies. The alternative expla- 
nation of the specific heat, that Ti is not a first-order 
transition, seems unlikely in light of not only the pres- 
ence, but also the shape of thermal hysteresis observed 
in p{T), which is particularly clear when a magnetic field 
is applied in the [110] direction (see Sec. IIVI below). 

The corresponding magnetic entropy 5, obtained by in- 
tegration, is shown in the inset of Fig.|Sl The two curves 
shown are calculated with two different extrapolations of 
Cp to zero temperature (dashed lines in the main panel). 
One extrapolation is obtained by connecting the Cp data 
point at the lowest T linearly to Cp = at zero tempera- 
ture, the other by connecting the Cp/T data at the low- 
est T linearly to Cp/T = at zero temperature. Given 
typically observed specific heats of Gd compounds, the 



actual specific heat at low temperature will most likely 
be between these two extrapolations, and thus the two 
curves shown in the inset of Fig. [S] may be considered 
lower and upper limits of the Gd 4/ entropy (neglecting 
additional experimental uncertainties). At 30 K the mag- 
netic entropy already reaches 80% (78% and 86% for the 
two low-T extrapolations made) of the full entropy of Gd 
4/ moments (R In 8, dashed line). Above 45 K it hardly 
varies anymore, having reached > 90% and > 97% of the 
full 4/ entropy for the two low-T extrapolations made. 
This strongly suggests that the Gd 4/ moments are not 
ordered above T2, and the transition at T2 corresponds 
to the ordering of the 4/ moments. 

By careful measurement and analysis of the low-field 
thermodynamic and transport data of solution grown sin- 
gle crystals of Gdl48 we have shown that there are sharp 
features associated with two magnetic phase transitions. 
These transitions appear to be primarily related to the 
ordering of the Gd sublattice. In addition we have been 
able to establish a clear correlation between small varia- 
tions in the Fe stoichiometry and T/v , Ti , and T2 . Having 
done this, we will now focus on the sample that is clos- 
est to full stoichiometry (crystal I) and try to more fully 
delineate and understand the field and temperature de- 
pendence and the structure of the magnetic phases in this 
sample. 



IV. HIGH-FIELD MEASUREMENTS 

Applying a magnetic field often helps to clarify the 
nature of magnetic phases observed in zero (or low) 
field measurements. We therefore measured the field 
dependent electrical resistivity and the magnetization 
on several samples. Shown in this section are electri- 
cal transport (subsection IIV A|) and magnetization (sub- 
section IIVB|) data measured on samples la/Ib at low 
temperatures in fields applied in-plane and out-of-plane. 
The features in resistivity corresponding to the two low- 
temperature transitions are much sharper when an in- 
plane field is applied and observed features suggest com- 
plex domain effects. The magnetization data are used 
to estimate the in-plane and out-of-plane components of 
the spontaneous ferromagnetic moment. The phase dia- 
grams for _ff II [110] and _ff || [001], which can be constructed 
from these transport and magnetization data, will be pre- 
sented and discussed in Sec. I VII 



A. Electrical transport 

Figure El shows the resistivity p (j||[110]) of sample la 
when a field of 3.5 T is applied either ||j (longitudinal re- 
sistivity) or II [110] (i.e. transverse resistivity, _Lj, but in 
a crystallographically equivalent in-plane direction). For 
both field directions, the magnitude of the jump at Ti is 
enhanced and the hysteresis in temperature greatly in- 
creased, making the first-order nature of the transition 



6 




-1 ' 1 ' 1 ' 1 ' 1 ' 1 

10 20 30 40 50 60 



r(K) 

FIG. 6: (Color online) resistivity p vs T in fields of 3.5 T ap- 
plied in-plane parallel or perpendicular to the current density 
J II [110] (sample la). Arrows indicate the T direction in which 
the measurements were performed. The open circle indicates 
the value of the Ti transition in 3.5 T, averaged for T increas- 
ing and decreasing. 



more apparent. The average (T jjT |) temperature of 
the Ti feature, however, is only weakly influenced by an 
in-plane field. In contrast to this, the T2 feature is shifted 
to higher T by -ff||[110], suggesting a ferromagnetic na- 
ture of the phase at T < T2. Furthermore, the fields 
remove the peak in p just above the drop at T2, and 
cause the feature to exhibit hysteresis. This different 
response to iJ||[110] confirms that there are indeed two 
independent magnetic transitions Ti and dividing the 
ordered state below T/v into three magnetic phases, a low- 
temperature phase (LTP), an intermediate-temperature 
phase (ITP), and a high-temperature phase (HTP). 

Whereas there are similarities in p{T) for the two field 
directions considered, there are also very pronounced dif- 
ferences, particularly for the Ti feature, below which the 
field cooled p{T) is much higher for iJ-Lj than for H\\j. 
A transverse magnetoresistance higher than the longi- 
tudinal one might arise from a bending of the electron 
trajectory by the Lorentz force (see, e.g., Ref. 123). How- 
ever, the zero-field-cooled p{T) below 33 K (black dotted 
line) is the same (within error bars) for H\\j and HL j. 
Indeed, as we will see (Fig. |S1 below), at 5K after cool- 
ing in zero field the resistivity is only very weakly field 
dependent for an in-plane field (both for transverse and 
longitudinal configurations) . 

This indicates that the difference in the p(T) curves of 
Fig. originates from the magnetism. The magnetore- 
sistance then may arise either from spin-disorder scat- 
tering or scattering of domain walls'^^ and the difference 
can be due to a different configuration of domains. In 
principle, this configuration (domain population, aver- 
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FIG. 7: (Color online) a) Resistivity p vs T (sample la, 
j||[110], T increasing) in fields applied parallel to ji||[110]. 
From top to bottom (thick line), 0.1, 0.2, 0.4, 0.5, 0.8, 1 
(thick line), 1.5, 2, 2.5, 3.5 T. b) As panel a), but in fields ap- 
plied parallel to [001] (j||[110]). From top to bottom (thick 
line), 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.8, 0.9, 1, 2, 5.5 T. 



age size of domains, etc.) may be different for i/||[110] 
(transverse magnetoresistance) and iJ||[110] (longitudi- 
nal magnetoresistance) due to different demagnetizing 
effects (the sample was cut in a rod such that the [110] 
direction is a bout three times longer than the [110] direc- 
tion). In any case, the fact that the resistivity is almost 
field-independent (for an in-plane field) at 5 K indicates 
that the relevant (for resistivity) domain configuration in 
the LTP is very hard to change by applying a magnetic 
field, whereas a field applied in the ITP prepares a cer- 
tain LTP domain configuration upon cooling below Ti. 
We note that application of a field in the ITP also seems 
to lead to small, but systematic, changes in the resistivity 
values in the HTP, persisting up to ^Tjy. 

Figure [7^) shows the evolution of p{T) [for clarity 
only the data measured upon warming, after the low- 
est temperature had been reached by cooling with the 
same field applied (field-cooled-warming measurement 




FIG. 8: (Color online) Resistivity vs ffUjH [110] at selected T 
(sample la). 

protocol) are displayed] with increasing H\\j (again sam- 
ple la, j II [110]). It seems clear that the T2 transition 
systematically shifts to higher temperatures with in- 
creasing H . The Ti transition is far less field depen- 
dent (for _ff||[110]), and as discussed above, quite hys- 
teretic. Furthermore, the resistivity in the ITP becomes 
lower in higher i?, developing into a well-defined weakly 
temperature-dependent valley above '^2 T. This develop- 
ment is also visible in the p{H) dependence at constant 
T (Fig.|HJ). There is only a minimal field dependence of p 
at 5 K in the LTP. In contrast to this, a pronounced drop 
in low fields at 25 K in the ITP is shifted to higher fields 
by increasing T and seems to transform into the drop 
associated with the ITP-HTP transition. In the p{T) 
curves of Fig. Et^), for intermediate field ranges, there 
are indications for an additional feature in the ITP-HTP 
transition. The details of this feature have to be clarified 
in future work. 

The evolution of p{T) with increasing i7||[001] for the 
same sample and current density direction is shown in 
Fig. [3)). The data shown were again obtained using 
the field-cooled-warming measurement protocol, but we 
note that application of H\\ [001] (unlike H\\ [110]) quickly 
suppresses any hysteresis, and broadens the transitions. 
Apart from a broadening, the T2 transition does not seem 
to be strongly affected. In contrast to this, i?||[001] sys- 
tematically shifts the Ti transition to lower temperatures 
and above ~1 T there are no longer any indications for 
a transition to the LTP. This indicates that the LTP is 
suppressed by -ff||[001]. 



B. Magnetization 

In order to investigate the magnetic nature of these 
phases, we measured the magnetization M as a function 
of T and H ([[[110] and ||[001]). In Fig. we display 
for sample lb M{H) loops at selected T. For i?||[110], 
at 40 K, M{H) shows almost linear behavior, but at 2K 
M{H) has a behavior typical for soft ferromagnets: satu- 
ration in a field which is of the order of the estimated de- 
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FIG. 9: (Color online) Magnetization Al as a function of field 
H at selected T, HjllllO] (blue) and //[[[OOl] (red), sample lb. 
Thin arrows indicate the H direction of the measurements, 
thick arrows indicate the criterion taken for the estimation of 
the spontaneous ferromagnetic moment. 



magnetizing field, but no hysteresis or remanence. How- 
ever, the saturation value of 4.5/iB/f.u. at 2K is much 
lower than the full Gd moment of 7pB /f .u. . At 26 K, pro- 
nounced hysteresis is visible in low H. This is ubiquitous 
for A4{H) in the ITP and correlates with the hysteresis 
in the drop in p visible in Fig. |H^) . The hysteresis may 
suggest a first-order metamagnetic transition. In higher 
fields M(H) is nearly constant, resembling the behavior 
at 2 K. In low fields M{H) indicates a spontaneous ferro- 
magnetic moment, unlike at 2 K, at 26 K we did observe 
remanence. The T dependence of the spontaneous ferro- 
magnetic moment, as estimated from Fig. taking into 
account demagnetizing effects (we took the magnetiza- 
tion at the first kink in low H, indicated by thick arrows 
in Fig.O, is shown in Fig. 1101 (black squares). The LTP- 
ITP transition in M{H), _ff||[110] is accompanied by a 
drop of the spontaneous moment to ^IpB^ and by the 
appearance of hysteresis. 

For -ff||[001] at 26 K (ITP), the form of M{H) (Fig. 
laired curves) is consistent with a spontaneous ferromag- 
netic moment of about 2 ps- For the LTP consider M{H) 
at 2K. Below 0.5 T, M{H) is linear, above it rises more 
strongly before slowly starting to saturate. The field Hi 
where M{H) starts to increase faster than linear is lower 
at higher T, extrapolating to zero at Ti. For samples 
with a higher Ti, Hi at low temperatures is proportion- 
ally higher (not shown). Whereas the feature at Hi as 
defined above is rather sharp, the slow saturation of M 
above Hi may be due to a small inhomogeneity of the Fe 
content in the sample, as discussed in Sec. IIIII 

For ff||[001], the saturation value of M < ^pb/^-Vl. 
is higher than for i?||[110], but still far below the full 
Gd 4/ moment. Note that the difference to 7/iB/f.u. is 
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FIG. 10: (Color online) Spontaneous ferromagnetic moment 
Ms, estimated from the M{H) curves in Fig.^lvs temperature 
T, i/||[110] (black squares) and //[[[OOl] (red circles), sample 
lb. The square root of the sum of the squares for both field 
directions is shown as open triangles. 



much too large to be explained by a negative contribu- 
tion of the Gd conduction electrons. The low saturated 
moment implies that the magnetization does not corre- 
spond to the magnetization of a simple ferromagnet. For 
i?||[001], the temperature dependence of the estimated 
spontaneous ferromagnetic moment is shown in Fig. 1101 
as red circles. It is similar to M{T) in low H (c.f. Fig. 

EJ. 

The linear M{H) behavior (below 0.5 T at 2K) may 
be due in-plane ferromagnetic moments resisting a reori- 
entation by an out-of-plane field. The LTP-ITP transi- 
tion could correspond to a reorientation of ferromagnet- 
ically aligned 4/ moments (as suggested in Ref. 0). Fig- 
ure ^1 indicates that in the LTP the ferromagnetically 
aligned moments are confined to the a—b plane, whereas 
in the ITP the moments are mainly out-of-plane, point- 
ing into a direction with an angle of about 25° to [001]. 
A crude estimate of the magnitude of the spontaneous 
moment is given by {M'^^-^q -\- M^qq;^)-^/^, shown in Fig. 
1101 as open triangles. Apart from one point near the 
LTP-ITP transition (where the estimation of Mg is less 
clear than in other regions) the temperature dependence 
of (M^^io + MloQ^y^^ is consistent with the LTP-ITP 
transition not affecting the magnitude of the ferromag- 
netic component, but only it's direction. Finally, we note 
that, like the low-field magnetization (Fig. ^ full line), 
the estimated spontaneous moment extrapolates to zero 
around 30 — 35 K, i.e., slightly higher than T2 as indi- 
cated by resistivity and the specific heat maximum. This 
roughly correlates with the onset of the 4/ contribution 
to the specific heat (Fig. (S)), which is not sharp at all. 
We therefore think it is not an indication of an additional 
transition, rather it may be due to an inhomogeneous Fe 
distribution or due to significant precursor-effects (e.g. 



fluctuations). 

The electrical-transport and magnetization data pre- 
sented in this section will allow us to construct H—T 
phase diagrams for -ff||[110] and -ff||[001]. We wiU do 
this in Sec. IVII Furthermore, we have established that 
the Gd 4/ moment order setting in at T2 has a large fer- 
romagnetic component, and the first-order (see Fig. O 
transition at Ti involves a reorientation of the ferromag- 
netic component into the a—b plane for T < Ti. In order 
to elucidate the underlying magnetic structure we will 
now focus on scattering data. 



V. X-RAY MAGNETIC RESONANT 
SCATTERING 

In order to prove the long-range nature of the magnetic 
order, and to learn more about the magnetic structure 
than the information that is obtainable from thermo- 
dynamic and transport measurements, scattering tech- 
niques can be applied, typically neutron scattering. For 
Gdl48, neutron scattering is less feasible because of the 
large neutron absorption cross-section of the natural iso- 
tope of Gd. An alternative technique, which relies on in- 
tense synchrotron radiation and exploits a large increase 
in the scattered intensity due to resonant processes when 
the energy of the radiation is close to an absorption edge 
of a magnetically ordered element, is x-ray resonant mag- 
netic scattering (XRMS)iSS*2Si2i To gain further insight 
into the magnetism of Gdl48, we therefore performed a 
first XRMS experiment on the (together with piece la 
from the same crystal) best characterized sample lb. 

Based on the neutron scattering and XRMS results on 
R148 with other rare earths,^i^i^i^°i" as well as on the 
thermodynamic and transport results presented in ear- 
lier sections, we expected at temperatures between T2 
and T^v a modulation of Fe moments with propagation 
along the [110] direction to be present, with a feasible 
resonant enhancement of the intensity scattered at the 
corresponding satellite reflections at the Gd Lu edge. 
Figure lllb ) shows reciprocal space scans along [h,h,0] 
with (7 — !■ TT polarization analysis and the energy tuned 
to the Gd Lu edge (7.934 keV) at 37 and 160 K. At 37 K, 
the charge reflection at (4,4,0) is accompanied by small 
satellites at (4 ± t,4 ± t,0) with r = 0.086. The satellite 
reflections are not present at 160 K, indicating that they 
are due to the magnetic ordering at the Neel tempera- 
ture (c.f. Fig. 121). We also found satellite reflections at 
(3 ± T,3 ± T,0), (2 ± r,2 ± t,0), and (1 ± r,l ± t,0), with 
0.06 < T < 0.135, depending on T. The fuU width at 
half maximum (FWHM) of the satellite reflections is ap- 
proximately twice the FWHM of the charge reflections. 
The corresponding smaller correlation length of the mag- 
netic order may be an indication of the presence of small 
domains. We roughly estimate the average size of the do- 
mains (in [110] direction) to be ~0.17/im or about 140 
unit cells. Rearrangements of these domains might be 
the cause of the small changes in resistivity values in the 
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FIG. 11: (Color online) a) Reciprocal space scans along [h,h,Q\ 
at 37 and 160 K (sample lb, Gd Lu edge). The data in the 
range h = 3.96 to 4.04 have been divided by a factor of 10000. 
b) Energy scans of the intensity of the (4 — r,4 — t,0) satellite 
reflections at 9, 37, and 130 K. Note that intensities in panel 
a) are not comparable to intensities in panel b). 



HTP (see the discussion of Fig. El in Sec. lIV^ . 

The intensity of the (4±T,4±r,0) satellite reflections is 
about four orders of magnitude lower than the intensity 
of the (4,4,0) charge reflection, despite the use of the po- 
larization analysis, which suppresses the charge reflection 
intensity by three orders of magnitude. Due to the small 
intensity of the satellite reflections and their vicinity to 
the charge reflections, there is considerable intensity due 
to the charge tail at the position of the satellite reflec- 
tions. A slight asymmetry in the charge tail leads to a 
larger charge contribution to the scattered intensity at 
the (4 -I- r,4 -I- r,0) reflection, and we therefore further 
analyzed only the (4 — t,4 — t,0) satellite reflection. 

FigurelTTb) shows, at selected T, energy scans over the 
Gd Lu edge for the satellite reflections at (4 — r, 4 — r, 0) . 
The intensity of the resonance is largest at 37 K. The 
lower intensity at 130 K is due to the vicinity of the Neel 
transition, the lower intensity at 9 K will be discussed 
below in Sec. IVfl The position and the (amplitude- 
normalized) peak shape of the resonance is T indepen- 
dent within our resolution and consistent with dipolar 
(£■1) transitionsiSSiSSiSi We note that a resonant en- 
hancement at the Gd Lu edge was observed only at the 
(4±T, 4±T, 0) satellite reflections, there was no enhance- 



ment at the (4,4,0) charge reflection, consistent with res- 
onant magnetic scattering. Reciprocal space scans and 
the resonant enhancement at the Gd Lu edge are consis- 
tent with a magnetic modulation with propagation vector 
(t,t,0) associated with Gd. 

The specific-heat measurement discussed in Sec. liiil 
implies that in the HTP the Gd 4/ moments are not 
ordered. Since XRMS is element-sensitive, this might 
seem in contradiction to our observation of a resonant 
magnetic signal at the Gd Lu edge at 37 and 130 K (Fig. 
\n\ . However, XRMS {El) of the Gd L edges is sen- 
sitive to the polarization of the Gd bd bands, not (di- 
rectly) to the 4/ moments iiS^iSS42I A Gd bd band polar- 
ization can also be induced by the Fe Zd moments. For 
this situation, a very weak intensity of the corresponding 
magnetic satellites is expectedi^ Compared to other Gd 
compoundsr2i2fli2ii2S the intensity of the satellite reflec- 
tions we observed is indeed very weak — more than three 
orders of magnitude lower. This low intensity might be 
due to the ordered moments pointing into a direction 
close to [001], since in our scattering geometry we are 
not sensitive to moments in [001] direction. However, 
changing the orientation of the sample such that [001] 
is in the scattering plane^*^ did not significantly change 
the intensity, suggesting that its low value is not related 
to the direction of magnetic moments. The behavior is 
similar to the situation in Dyl48 at high temperaturesii 

At this point it may be useful to briefly comment 
on an XRMS signal coming directly from the Fe mo- 
ments. This would of course be very important to con- 
firm the above considerations. Unfortunately, the en- 
ergies of the L absorption edges of id transition met- 
als do not allow diffraction experiments for regular crys- 
talline materials Resonances at the K edge may be 
used instead, but the x-ray resonant process at transi- 
tion metal K edges is neither well understood nor ef- 
ficient. Reported resonant enhancements at transition 
metal K edges of magnetic scattering intensities were 
typically v ery small and often broad (see, e.g., Refs. 
13411851186118711881 ) . compared to the substantial enhance- 
ments often observed at rare earth L edges (see for ex- 
ample Refs. I29l3nl3ll32tl . 

An additional problem in our case is the large charge 
scattering background due to the vicinity of the magnetic 
satellite reflections to the charge reflections. This is more 
problematic for Fe K edge measurements, because at the 
Fe K edge, the angle of the analyzer crystal for the polar- 
ization analysis is further away from the ideal 90° than 
in the case of the Gd Lu edge, resulting in a significantly 
higher charge scattering background at the Fe K edge. 
Performing energy scans in the range of the Fe K edge 
(7.11 keV) we did not observe a clear resonant enhance- 
ment. In view of the mentioned difficulties, this does not, 
however, imply that the Fe Zd moments are not ordered. 
Rather, we conclude that the direct observation of an 
XRMS signal from the Fe moments will need a separate 
experiment with specialized equipment and enough time 
to gather high statistics. 
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Although we did not observe a clear resonant enhance- 
ment at the Fe K edge, based on the above discussion 
of the intensity of the XRMS at the Gd Ln edge and 
the comparison with other R148 compounds, we conclude 
that the (r, t, 0) propagation vector describes the order 
of the Fe moments. A Fe moment modulation along [110] 
is common for R148 compounds, but in R148 with other 
rare earths r is typically larger (0.127 < r < 0.243) and 
varies less with T.^i^i^°i" 

At 9 K, the satellite reflections are present as well [Fig. 
Illb )]. with the propagation vector in the same direc- 
tion. This demonstrates that the Fe moments maintain 
their modulation from the HTP, where the Gd 4/ mo- 
ments are not ordered, in the LTP, where the Gd 4/ 
moments are ordered. This is similar to the situation 
in other R148 compounds. However, the low intensity 
of the (4 — T, 4 — r, 0) satellite reflections at low T indi- 
cate that for Gdl48, unlike the situation in other R148 
compounds, the rare-earth 4/ moments do not fol- 
low the Fe moment modulation. In order to find possible 
modulations involving the Gd 4/ moments, we searched 
at 9 K also for magnetic satellite reflections in other high- 
symmetry directions. However, full scans along [110], 
[100], and [001] (for scattering geometries mentioned 
above) did not turn up any reflections with intensity high 
enough to correspond to a modulation of a significant 
part of the 4/ moments. 

Further XRMS experiments, including tracking the 
observed magnetic satellite reflections in detail through 
temperature and applied field, and comparing the res- 
onances at absorption edges other than Gd L\i are 
planned. 



VI. DISCUSSION AND CONCLUSIONS 

In the following we will discuss our results, starting 
with the magnetic ordering at the Neel temperature and 
the nature of the resulting phase. We will then discuss 
the H — T phase diagrams (for two field directions) at 
low temperature and the magnetic order in these phases. 

Our resistivity and magnetization data indicate that 
GdFe4Al8 orders magnetically at « 155 K (see Sec. 
I1II|I . Even small (1 — 3%) iron deficiency significantly 
lowers T/v, indicating that magnetic interactions involv- 
ing the Fe 3d moments drive the transition. In the x-ray 
resonant magnetic scattering (XRMS) experiment (Sec. 
IVj) . the order below Tn was identified as an incommen- 
surate antiferromagnetic order with propagation in [110] 
direction. Although corresponding satellite reflections 
were measured on the Gd Lu edge, the T dependence of 
entropy associated with the Gd 4/ moments, estimated 
from specific-heat data (see Fig. |S1 in Sec. IIII|I , indicates 
that at high temperature, the Gd 4/ moments do not 
participate in this antiferromagnetic order, which is thus 
driven by Fe-Fe moment interactions. This can explain 
the weak intensity of the measured satellite reflections. 
Antiferromagnetic order propagating in [110] direction. 
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FIG. 12: (Color online) Low-temperature phase diagrams of 
GdFe3.96(i)Al8 (sample la/b) for -ff||[110] (a) and -ff||[001] 
(b). Triangles in a) are from resistivity measurements on 
sample la, the other symbols from magnetization measure- 
ments on sample lb (see text for details) . The three magnetic 
phases LTP (low-temperature phase), ITP (intermediate- 
temperature phase), and HTP (high-temperature phase), as 
well as the zero-field transition temperatures Ti and T2 are 
indicated in both panels. A possible fourth magnetic phase 
existing in fields applied parallel to [110] (HFP) is indicated 
in panel a) . Underlaid is a contour plot of the magnetization 
M (sample lb) . 



at high T confined to Fe 3d moments, seems to be rather 
typical for i?Fe4Alg compounds. 

In addition to the Neel transition, wc found two more 
phase transitions at low temperature (Ti and T2), di- 
viding the ordered state into three phases: LTP (low- 
temperature phase), ITP (intermediate-temperature 
phase), and HTP (high-temperature phase). Due to the 
proximity of the two low-T transitions it proved to be 
very useful to apply magnetic fields in two directions and 
thereby more clearly separate these phases. 

Low-temperature phase diagrams established from re- 
sistivity and and magnetization measurements (see Sec. 
IIV|) on crystal I (samples la and lb) are presented in 
Fig. US [for i?||[110] in panel a) and for i?||[001] in panel 
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b)]. The phase diagrams are underlaid with a contour- 
plot of the magnetization [data determined from M{H) 
measurements as displayed in Fig. El for H increasing] . 

For -ff||[110] [Fig. nil)] the phase boundaries (A,v) 
were determined from p{T) measurements [steepest 
slope, c.f. Fig. [7^); shown is the average between points 
determined on curves measured with T increasing and 
decreasing]. In M{T,H) the ITP-HTP transition corre- 
sponds to, upon decreasing T, a steep rise in M and satu- 
ration. For definiteness, we used the extrema in AM/dH 
[d(Mr)/dT], plotted in Fig. HSli) as o (□). In case of 
hysteresis we again took the average between T or H 
increasing and decreasing. Above 0.5 T, this tracks the 
ITP-HTP transition deduced from p(T) well, but in lower 
H the smallness of the magnetization component in [110] 
direction leads the transition line established by these cri- 
teria two lower T, in fact extrapolating to Ti for H = 0. 
The M{H) curve measured at 26 K (Fig.lHI) and the p{H) 
curve measured at 25 K (Fig. IS)) do suggest an additional 
metamagnetic transition in the temperature region of the 
ITP. This transition is what is determined from the crite- 
ria above below 0.5 T, starting from Ti in zero field and 
apparently merging with the boundary to the HTP at 
about 28 K in 0.5 T. 

This seems to separate the ITP as deduced from resis- 
tivity measurements into two phases, we call the tentative 
high field phase HFP. In the {longitudinal, iJ||[110]) mag- 
netization the HFP is indistinguishable from the LTP, 
and thus it could be that the HFP is better described 
as a "modified LTP" phase rather than as a "modified 
ITP" one. However, the electrical transport data (see, 
e.g., Fig.EJ clearly indicate a first-order transition sepa- 
rating LTP and HFP. This transition exhibits a wide hys- 
teresis in temperature when i/||[110] is substantial (see 
Fig. El and likely involves complex transformation of do- 
mains present in the LTP and in the HFP (see the brief 
discussion in Sec. IIV Al note that the HFP was not yet 
introduced there) warranting further investigation. The 
connection of the LTP-HFP transition with pronounced 
domain effects might indicate that domains are a signif- 
icant force driving this transition, but presently this is 
tentative at best. The large hysteresis associated with 
the LTP-HFP transition complicates the determination 
of the equilibrium transition line. Estimating the equilib- 
rium transition by taking the average of the transitions 
for T increasing and decreasing [for 3.5 T this is indicated 
by an open circle in Fig.|Hl in the phase diagram, the tran- 
sition is indicated by open down triangles] suggests that 
i/||[110] favors neither the LTP nor the HFP strongly 
Similarly, the transition between the ITP and the HTP is 
hard to pin down from M{T) measurements with _ff || [110] 
because of the smallness of the component along [110] of 
the magnetization. We observed no second "extrcmum 
in [d(MT)/dT]", using this criterion only leads to the 
transition between the HFP and the ITP, rather than 
to the one between the ITP and the HTP. For a field 
applied along [110], resistivity vs temperature measure- 
ments clearly indicate two transitions, whereas (longitu- 



dinal) magnetization vs temperature measurements in- 
dicate only one. Only when examining the spontaneous 
magnetization along [110] fFig. llOl) is the presence of two 
separate transitions clearly indicated. 

For i?||[001] [Fig. [13))], three distinct regions in the 
contour-plot of the magnetization are readily identified 
with LTP (low M), ITP (high M) and HTP (low M). 
The transition lines between these phases arc qualita- 
tively consistent with what may be deduced from the re- 
sistivity measurements presented in Fig.Et>)- In contrast 
to the phase diagram for _ff||[110], there are no indica- 
tions for the presence of additional phases for 7J||[001]. 
Since H\\ [001] significantly broadens the features in p(T), 
we rely here on magnetization measurements (see Figs. 
21 and O only to draw the transition lines in Fig. 112b'). 
To fix the ITP-HTP transition, we used the same proce- 
dure as for ffll[110] (squares). For i/||[001] too, the H 
dependence of the transition and a saturation consistent 
(considering demagnetizing effects) with a spontaneous 
moment indicate a ferromagnetic component in the ITP 
(see also Fig. lTUI) . The LTP-ITP transition is determined 
from the M{H) loops (c.f. Fig.O, taking the field value 
where M starts to increase faster than linear as the crite- 
rion for the transition field. This transition is displayed 
by yellow circles in Fig. 112b). For H ^ the transition 
line extrapolates to Ti, reaffirming it's identity as the 
LTP-ITP transition. It can be seen that i?ll[001] quickly 
suppresses the LTP. 

Based on the specific-heat data, we can conclude that 
the HTP to ITP transition corresponds to the ordering of 
the Gd 4/ moments. Since the magnetization data indi- 
cate the presence of a spontaneous ferromagnetic moment 
in the ITP, this Gd 4/ ordering is cither ferromagnetic 
or has a ferromagnetic component. The increase of the 
transition temperature T2 upon application of a magnetic 
field fFig.ll2|) is consistent with this. Thus, Gd moments 
do not simply follow the Fe moment order from higher 
temperatures [this is in contrast to, e.g., DyFe4Al8 (Ref. 

a. 

The transition from the LTP to the ITP seems to be 
of first-order as indicated by the observation of thermal 
hysteresis in the resistivity (see Fig. O as well as the 
magnetization (see Fig.j^J. Surprisingly though, no fea- 
ture could be identified in the specific heat corresponding 
to this transition and any latent heat would have to be 
very small. The LTP-ITP transition is associated with a 
lock-in of the observed spontaneous ferromagnetic com- 
ponent fFig. llOfl into the a—b plane in the LTP, but there 
are no clear indications that its magnitude is changed by 
the transition. The appearance of an additional antifer- 
romagnetic modulation is possible, at least it is not clear 
why a pure lock-in transition of a ferromagnetic compo- 
nent should result in a sizeable increase of the resistivity 
(both for j||[110] and jl|[001]). 

The presence of the HFP in i?||[110] complicates the 
situation. In fact, there are clear indications (meta- 
stability) for the first-order nature of the transitions at 
all boundaries of the HFP. In contrast to this, there are 
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no indications for meta-stability in the ITP-HTP transi- 
tion. Meta-stability seems to be present in the LTP-ITP 
transition, but it is much weaker than in any transition 
from or to the HFP. We point out that the first-order 
transition fines LTP-HFP and HFP-HTP both seem to 
extrapolate to Ti in zero field, and therefore the observed 
smaller meta-stability of the LTP-ITP transition might 
well be connected with the HFP phase. This is a strong 
incentive for further investigations of the HFP. 

We now attempt to draw a picture accounting for the 
full magnetic moment of both Gd and Fe. The full or- 
dered moment of Gd is T/is/Gd. Since the valence of 
Gd is fixed and there are no crystal-electric-field effects, 
this is the same for different compounds, with only small 
corrections, up to ~ 0.5 /is, due to contributions by itin- 
erant electrons (such as the Gd 5d electrons). The full 
moment of Fe can vary from compound to compound. 

We know from the specific heat and XRMS measure- 
ments that in the HTP Gd moments (the Gd 4/ mo- 
ments at least) are not ordered, and Fe moments are an- 
tiferromagnetically ordered with propagation along [110]. 
Magnetization measurements further show that there is 
no significant ferromagnetic component present [see the 
40 K M{H) loops in Fig. We recall that, according 
to the x-ray-diffraction-structure refinement and specific- 
heat measurements, the Gd lattice is fully occupied and 
that at low temperatures the Gd 4/ moments are well or- 
dered. The appearance of a spontaneous ferromagnetic 
moment below the temperature T2 where the Gd mo- 
ments start to order (compare Figs. 13 and indicates 
that the ordering of the Gd moments is ferromagnetic, 
but the saturation moments for both field directions are 
far below the T/is/f.u. (see Fig. O expected for a full 
ferromagnetic alignment of Gd moments. What could be 
the reason for this? 

One possibility is that the Gd 4/ moments are not 
fully ferromagnetically ordered, but rather also have an 
antiferromagnetic component. However, since we did not 
find any satellite reflections with high enough intensity 
to account for the necessary fraction of the full Gd 4/ 
moment in any of the high symmetry directions (see Sec. 
IVl note that the intensity of the (4 — r, 4 — r, 0) reflection 
is much too low) this scenario is not very likely. 

From high-field measurements on finely ground free 
powder particles of Gdl48 (likely, at least roughly, cor- 
responding to iJ||[001], ITP), Duong et aim- suggested 
for iJ — » a ferrimagnetic structure with Gd moments 
parallel, Fe moments antiparallel to the field, with an an- 
tiferromagnetic Fe moment order re-established only in 
high fields. On the one hand, such a simple ferrimag- 
netic arrangement of Gd and Fe moments is inconsistent 
with our observation in the XRMS experiment of the 
presence of the (4 — r, 4 — r, 0) satellite refiection charac- 
teristic for the antiferromagnetic Fe modulation even at 
9K [Fig. II lb)]. On the other hand, a related idea would 
be consistent with both the XRMS and the magnetiza- 
tion measurements: upon the ferromagnetic ordering of 
the Gd 4/ moments, the Fe moments could pick up a 



ferromagnetic component, leading to a canted antiferro- 
magnetic structure. Such a canted magnetic structure 
with a large ferromagnetic component may also explain 
why the intensity of the (4 — r, 4 — r, 0) satellite refiection 
is about three times smaller at 9 K (in the LTP) than at 
37 K (in the HTP). Furthermore, it would also allow Fe 
to influence the low-T transitions, as indicated by the 
influence of Fe deflciencies on all characteristic temper- 
atures (see Sec. IHI|I . The results of a recent dichroism 
experiment (on a crystal ground to powder), which will 
be discussed in detail elsewhere, are consistent with such 
a scenario as well. 

This picture covers, at least for the LTP and the HTP, 
all the data presented in the paper, but without more de- 
tailed scattering data it remains speculative. Particularly 
concerning the ITP and HFP, the absence of scattering 
data makes detailed conclusions elusive. 

Since the direction of the spontaneous ferromagnetic 
moment is different for the ITP and the LTP, magne- 
tocrystalline anisotropy (MCA) is likely one of the forces 
driving this transition. In rare earth compounds, MCA 
typically arises from crystal-electric-field effects. As we 
mentioned in Sec. 12 this source of anisotropy is absent 
for Gd. Signiflcant MCA that was nevertheless observed 
for both elemental Gd and Gd compounds has been 
ascribed to a combination of dipole-dipole interactions 
and (mainly) anisotropy in the 5d bandsi224£ The corre- 
sponding magnetic-anisotropy energies are much lower 
than the magnetic-anisotropy energies resulting from 
crystal-electric-field effects in other rare earths. A weak 
magnetic-anisotropy energy might also, at least partly, 
explain why no latent heat could be observed at Ti. In 
Gdl48, anisotropy in the Fe 3d bands might be an addi- 
tional source of MCA4i In any case, the different spon- 
taneous moment direction in the ITP likely results in a 
different structure of the Fe moments. The apparent shift 
of Ti upon changing the Fe stoichiometry may be taken 
as an indication of the involvement of Fe moments in 
the LTP-ITP transition. Before going any further in the 
discussion of the driving forces (there is likely a delicate 
balance between multiple energy scales) of the transitions 
between LTP, ITP, and HFP, the magnetic structures of 
these phases need to be solved. With this aim in mind, 
additional scattering experiments are planned. 

Regardless of the details of the magnetic structure at 
low temperature, it is important to point out that i) on 
the one hand, the Fe moments have at least a compo- 
nent antiferromagnetically modulated in [110] direction 
ii) The Gd 4/ moments, on the other hand can not have a 
component ordered with the same propagation as the Fe 
3d moments, because then the intensity of the reflections 
measured at 9 K would have to be orders of magnitude 
higher (see Sec. 01. This is strikingly different from the 
behavior for example in Dyl48, where the rare earth mo- 
ments order at lower temperature, but then follow the Fe 
moment modulation, and it implies a co- existence of two 
distinct orders associated with Gd 4/ and Fe 3d moments. 
Whereas compounds in which one magnetic sublattice 
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orders at temperatures much lower than the other mag- 
netic sublattice are not very rare (indeed this was well 
known to be the case for R148 compounds in general), 
we are aware of only one example, PrBa2Cu3 06.92,'^^ in 
which the magnetic sublattice ordering at lower T has 
no component with the same modulation that the other 
sublattice has at all T. 

The at low temperatures persisting antiferromagnetic 
modulation of the Fe moments coexisting with an order of 
the Gd moments with a large ferromagnetic component 
and without a component having the Fe moment modu- 
lation is the central point of our paper. It may indicate 
a modified relation between the strength of the Gd-Gd 
and Gd-Fe moment interactions in Gdl48 as compared 
to other R148 compounds. As we mentioned, apart from 
the stronger couplings between 4/ and conduction elec- 
trons in general, the most peculiar property of Gd is the 
absence of crystal-electric-field effects. 

It seems remarkable that in the absence of the MCA 
due to crystal electric field effects, the magnetic behav- 
ior in Gdl48 is more complex than in other R148. Likely, 
strong crystal-electric-field effects play a central role in 
the magnetism in R148 (except for Gdl48). E.g., in 
Dyl48, crystal-electric-field effects press the 4/ moments 
into the crystallographic plane, making an order follow- 
ing the Fe-moment modulation more probable, since the 
Fe moments are already in-plane. Our results on Gdl48, 
in contrast, indicate a complex interplay of local 4/ and 
itinerant Zd moments when crystal-electric-field effects 
are absent. 



VII. SUMMARY 

We presented an extensive set of data including magne- 
tization, electrical transport, specific heat, and x-ray res- 
onant magnetic scattering, measured on flux-grown sin- 
gle crystals of GdFe4Al8. We found that two transitions 
at Ti ~ 21K and r2 ~ 27K at low temperature divide 



the ordered state below T/v 155 K, where Fe moments 
order antiferromagnetically with modulation along [110], 
into three phases. The corresponding phase diagrams for 
fields applied in [110] and [001] direction were presented. 

Gd 4/ moments order, mainly ferromagnetically, be- 
low T2. Above Ti the ferromagnetic component of the 
moments points into a direction close to [001], whereas 
below Ti the ferromagnetic component is locked into the 
a—b plane. The Gd 4/ moment order is distinct from 
the Fe-moment order At low temperature T: the Fe mo- 
ments still have at least a component modulated along 
[110], but the Gd 4/ moments apparently do not have 
such a component. We proposed a canted antiferromag- 
netic structure of the Fe moments at low T, which can 
cover all the low-T data presented in the paper, but to 
some extent still is speculative. Above The transition at 
Ti likely involves a delicate balance of multiple energy 
scales associated with both Gd and Fe. 

The complex magnetism in GdFe4Al8 as compared to 
i?Fe4Alg compounds with other rare earths R likely is 
related to a modified ratio of coupling strengths and to 
the absence of crystal-electric-field effects. In order to 
further elucidate the structure of the magnetic phases of 
GdFe4Al8 additional scattering experiments are planned. 
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